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Abstract 


In  this  report,  short  cylindrical  shells  are  studied 
under  pure  bending.  An  edge  kinematic  constraint  loading 
method  is  adapted  for  use  with  the  STAGSC  computer  code. 

With  it,  short  unstiffened  cylindrical  shells  are  studied 
for  prebuckling  displacements  and  compared  with  axially 
stressed  cylinders.  End  rings  are  added  to  nc  ^  their  in¬ 
fluence.  The  bending  is  applied  to  short  internally  stiff¬ 
ened  (smeared  stringer  theory)  cylinders,  with  the  result¬ 
ing  buckling  values  and  patterns  contrasted  against  axial 
loading.  The  L/R  range  of  interest  is  from  .698  to  1.361. 
Symmetrical  and  unsymmetrical  12  inch  cuts  are  added  to  the 
stiffened  cylinder  to  study  the  resulting  effects  from  this 
imperfection.  The  results  indicate: 

a.  For  short  length  (L/R  <  10)  cylinders,  end  loading 
by  pure  bending  causes  structural  deformation  patterns 
different  than  long  cylinder  bending  deflections. 

b.  The  imposed  boundary  effect  for  unstiffened  cylin¬ 
ders  diminishes  with  L/R  increase  until  L/R  =  20  where  there 
is  no  noticeable  deformation  influence. 

c.  End  ring  stiffeners  studied  did  not  alter  the  basic 
short  cylinder  deformation  patterns. 

d.  Compressive  zone  buckling  failure  under  bending  can 
be  modeled  through  axial  compression  analysis. 
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A  STUDY  OF  SHORT  CYLINDRICAL  SHELLS  WITH  AND 


WITHOUT  CUTOUTS  UNDER  PURE  BENDING 


I.  Introduction 


Prologue 

Anyone  having  bent  a  can  or  slender  tube  has  wit¬ 
nessed  cylindrical  shell  instability.  More  important  mani¬ 
festations  include  kinking  of  ocean  floor  piping,  nuclear 
plant  fluid  pipe  failures,  and  the  failure  of  aircraft 
fuselages  or  spacecraft  boosters.  The  aerospace  design 
engineer  is  faced  with  trying  to  optimize  the  strength  to 
weight  characteristics  of  these  flight  structures.  This 
compels  him  to  understand  the  problems  of  thin  shell  sta¬ 
bility  and  methods  for  determin'ng  design  characteristics. 
It  was  this  sort  of  impetus  that  motivated  Donnell  (Ref  16) 
to  develop  his  simplified  equations  of  equilibrium  for  a 
cylindrical  wall  in  1933.  The  high  speed  computer  has 
greatly  enhanced  the  ability  to  analyze  and  understand  thin 
shell  structures.  I  will  use  one  such  computer  code 
(STAGSC)  to  analyze  the  effects  of  pure  bending  on  short 
cylinders. 

The  thesis  topic  arose  from  discussions  with 
Dr.  Palazotto.  A  previous  student,  Norton  Compton  (Ref  15) 
had  problems  in  keeping  the  end  plane  from  warping  when  he 
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attempted  to  apply  pure  bending  moments.  My  task  became 
one  of  finding  a  method  for  applying  pure  end  rotation 
while  keeping  the  end  plane  undistorted  and  circular.  Once 
the  method  was  perfected,  I  was  to  study  its  effects  on 
stiffened  and  unstiffened  cylindrical  shells  with  and  with¬ 
out  cutouts.  Through  an  associate.  Dr.  Palazotto  received 
a  method  devised  by  Jorgen  Skogh  and  Frank  Brogan  (Ref  29) , 
which  I  adapted  for  STAGSC  (Ref  5) .  This  paper  will  cover 
previous  work  in  bending  analyses,  the  Skogh  routine,  and 
results  of  employing  it. 

Background  Theory  and 
Literature  Search 

A  shell  structure  is  usually  defined  as  a  body 
enclosed  between  two  closely  spaced  curved  surfaces. 
Cylinders  are  a  subset  in  that  there  is  only  one  radius  of 
curvature.  Shell  theory  is  an  attempt  to  use  the  assump¬ 
tion  of  smallness  in  the  thickness  direction  to  relate  the 
three-dimensional  problem  of  deformation  of  the  body  to 
that  of  displacement  components  on  a  middle  surface  or 
reference  surface.  The  analysis  therefore  involves  making 
assumptions  about  the  shells'  dimensions.  For  thin 
cylindrical  shells  the  following  assumptions  given  by  Love 
are  standard  (Ref  25) : 

1.  The  shell  is  thin.  This  means  that  the  thick¬ 
ness  of  the  shell  is  small  compared  with  the  radius  of 
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curvature  R  of  the  middle  surface,  so  that  the  ratio  is 
small  compared  to  unity. 

2.  The  deflections  of  the  shell  are  small,  and  the 
strains  in  the  direction  of  the  normal  are  small  enough  to 
be  neglected.  This  allows  us  to  refer  the  analysis  to  the 
initial  configuration  of  the  shell. 

3.  The  normal  stresses  acting  on  planes  parallel 
to  the  middle  surface  are  negligible  compared  with  other 
stress  components  and  may  be  neglected  in  the  stress-strain 
relations.  This  assumption  will  generally  be  valid  except 
in  the  vicinity  of  highly  concentrated  loads. 

4.  The  components  of  the  displacements  are 
linearly  distributed  across  the  thickness. 

5.  The  shear  strains  which  cause  the  distortions 
of  the  normals  to  the  middle  surface  can  be  neglected. 

Which  says  that  normals  to  the  undeformed  middle  surface 
remain  normal  to  it  after  deformation. 

Assumptions  2,  3  and  5  are  the  shell  counterparts  for  the 
Kirchhoff  assumptions  of  thin  plate  theory  (Ref  14).  For 
a  complete  listing  of  the  assumptions  for  Kirchhoff' s 
theory  of  plates  and  Love's  first  approximation  to  the 
theory  of  shells,  see  Gould  (Ref  17).  These  assumptions 
are  applied  to  a  portion  of  a  shell  surface.  Fig  1,  along 
with  the  internal  forces  and  moments,  expressed  as  forces 
and  moments  per  unit  distance  (intensities)  along  the  edge. 
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Figure  1.  Cylindrical  Shell  Segment 


On  Fiq  1,  the  and  Ng  are  in  plane  normal  intensities; 

NxQ  and  NQxare  shear  intensities;  Qx  and  QQ  are  transverse 
shearing  force  intensities;  the  Mx  and  M0  are  bending  moment 
intensities;  and  the  Mx0  and  M0x  are  twisting  moment  inten¬ 
sities.  See  Appendix  A  for  their  definitions  along  with 
the  constitutive  equations  and  kinematic  relations.  Allow¬ 
ing  a  slight  deformed  position,  the  forces  are  summed  and 
recombined  with  the  constitutive  and  middle  surface  kine¬ 
matic  relations  to  get  equilibrium  eqs  with  dead  load  pres¬ 
sure  loading  p: 


RNx,x  +  Nxe,e=  0 


dv*„  +  i  Ng 


RNxe,x  +  Ne,e  *  0  (1> 

<Nx”'xx  +  I  Nxow'xe+RiIN0w'ee>  -  P 
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where  V4w  =  w. 
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and  D  is  the  bending  stiffness  D  = 


Et 


12  (1-v2) 


Equations  (1)  are  coupled  nonlinear  partial  differential 
equations  for  quasi-shallow  cylindrical  shells.  Linear 
equations  are  obtained  by  deleting  all  quadratic  and  higher 
order  terms  in  u,  v,  and  w  since  they  are  products  of 
already  assumed  small  quantities.  The  resulting  equations 
are : 


RN  +  N  .  .  =  0 
x,x  x8,0 


RNX0,X  +  N0,0  =  0 


DVV  +  ±  N0 


=  p 


(2) 


The  difference  between  Eqs  2  and  those  for  a  flat  plate  is 
that  Eq  2c  is  not  uncoupled  from  Eq  2a  and  2b.  This  gives 
us  the  important  result  that  axial  and/or  circumferential 
movement  produces  radial  displacements. 

Eq  1  can  be  derived  through  the  use  of  the  minimum 
Potential  Energy  Criterion.  A  review  of  energy  methods 
and  variational  calculus  can  be  found  in  Brush  (Ref  14)  and 
Saada  (Ref  25) .  The  total  potential  energy  V  can  be 
denoted  as  the  sum  of  the  strain  energy  U  of  the  cylindri¬ 
cal  shell  and  potential  energy  0  of  the  applied  forces. 

The  condition  for  equilibrium  is  arrived  at  by  taking  the 
first  variation  of  the  total  potential  energy  6V  and  setting 


it  equal  to  zero,  while  the  criterion  for  stability  is  that, 
in  addition,  the  second  variation  62V/2  must  be  positive 
definite.  For  a  cylindrical  shell  the  strain  energy  U 
can  be  seen  as  a  combination  of  a  membrane  strain  energy 
Um  and  bending  strain  energy  U^. 

um  =  ¥  //{e£  +  zl  +  2vV:e  +  ^  Yie}  dxde  (3) 

Ub  =  ¥  //{kx  +  *6  +  2vicxKe  +  2(1~v)k:x0}  dxde  (4) 

Et 

where  C  is  the  extensional  stiffness  C  =  ^  v2j  • 

And  for  the  shell  subjected  to  dead-load  lateral  pressure  p, 
the  potential  energy  is  given  by 

=  -R//pwdxd0  (5) 

Setting  the  first  variation  6V  equal  to  zero  leads  again 
to  Eqs  1. 

Before  determining  the  buckling  equations,  a  quick 
overview  of  stability  theory  is  appropriate.  Instability  is 
a  nonlinear  phenomenon  in  which  a  small  increase  in  load 
causes  a  disproportionately  large  increase  in  deformation. 
This  load  is  called  the  critical  load  Pcr«  The  nonlinearity 
can  be  caused  by  either  material  (plasticity)  or  geometric 
nonlinearity.  Only  geometric  nonlinearity  is  considered 
in  this  report.  Sobel  (Ref  30)  details  how  the  loss  of 
stability  can  happen  through  either  bifurcation  buckling  or 
geometric  collapse.  Figure  2a  represents  bifurcation 
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buckling.  The  load  displacement  curve  emanating  from  the 
origin  is  intersected  by  another  load  displacement  curve 
that  represents  a  different  equilibrum  configuration.  The 
point  of  intersection  is  called  the  bifurcation  point,  and 
the  load  at  that  point  is  called  the  bifurcation  load. 
Figure  2b  represents  collapse.  The  load  displacement  is 
nonlinear  from  the  origin  and  reaches  a  relative  maximum. 
The  load  at  this  point  is  called  the  collapse  load. 


Lateral  Displacement  Lateral  Displacement 

Figure  2a.  Bifurcation  Figure  2b.  Collapse 

The  adjacent-equilibrium  path  shown  in  Fig  2a  is 
sloping  upward  indicating  a  stable  post  buckling  path. 

This  is  not  always  so  and  in  most  cases  does  not  happen. 
Koiter  (Ref  18)  studied  the  post  buckling  equilibrium  con¬ 
figurations  through  the  energy  formulation  of  the  loaded 
structure.  The  potential  energy  expression  used  earlier 
is  of  the  general  form 

V  =  /A  F(u,v,w,X)dA  (6) 
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where  u,  v,  and  w  are  middle  surface  displacement  components 
and  X  is  an  applied  load  parameter. 

Let  the  variables  become 

u  ■>  u  +  u 
o  1 

v  •+  v  +  v  (7) 

o  1 

W  -*•  w  +  w 
0  1 

where  (u  ,  v  ,  w  )  represents  an  equilibrium  position  on  the 
0  0  0 

primary  path,  (u,  v,  w)  form  a  nearby  adjacent  equilibrum 
position  corresponding  to  the  same  value  of  X ,  and  (u^ ,  v  , 
w^)  is  a  small  inf initisimal  displacement.  Then  when  the 
change  in  potential  energy  is  written 

AV  =  ^62V  +  i<$3V  +  (8) 

The  first  order  term  drops  out  since  (uo,vo,wQ)  is  an  equi¬ 
librium  position.  Application  of  the  stationary  potential 
energy  criterion  to  the  expression  Av  leads  to  the  equations 
governing  post  buckling  behavior.  The  cases  shown  in  Fig  3 
arise  depending  on  the  characteristics  of  the  structure. 

As  is  seen  in  Fig  3  the  post  buckling  paths  vary  consider¬ 
ably.  The  solid  line  secondary  path  represents  perfect 
structure  adjacent  equilibrium,  whereas  the  dash  dot  line 
denotes  imperfection  sensitivity.  Notice  that  the  imperfec¬ 
tion  sensitive  lines  are  nonlinear  from  the  origin  and  hence 
more  closely  model  collapse  modes.  In  all  cases,  the  imper¬ 
fections  lower  the  load  at  which  the  structure  becomes 


Case  I  Case  II  Case  III 

Stable-Symmetric  Unstable-Symmetric  Asymmetric 
( a '  ( b )  ( c ) 

Figure  3.  Load  Displacement  Curves  For  Perfect  and 

Imperfection  Sensitive  Structures.  From(^o) 

unstable.  It  is  most  critical  in  case  II  where  the 
imperfect  structure  will  snap  at  the  limit  point. 

The  equilibrium  equations  are  now  transformed  into 
a  linear  coupled  set  representing  the  buckling  mode.  First, 
the  reference  surface  displacement  components  are  replaced 
by  Eq  7,  in  Eq  1.  Accordingly,  the  force  intensities 
become 


N 

-+  N  + 

AN 

X 

xo 

X 

Ne 

-  N0o  + 

AN0 

(9) 

N 

X 

"  NX0O+ 

*Nxe 

Where  terms  with  o  subscripts  refer  to  the  primary  equili¬ 
brium  configuration,  and  ANx,  ANg,  and  ANx0  correspond  to 
the  change  due  to  increment.  When  Eqs  7  and  Eqs  9  are 

substituted  into  eqs  1,  all  terms  in  u  ,  v  ,  w  alone  as 

o  o  o 
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well  as  terms  such  as  Nxq  fall  out.  Also,  if  the  higher 
order  terms  in  u,  v,  w  and  etc.  are  neglected  due  to 
their  small  relative  size,  the  remaining  equations  become 


RN  .  "^NQ1Q  =  0 
xl,x  X01,0 


RNX01,X  +  N01,0  =0 


(10) 


dV-w.  +  .  -  |(N  w.  +  w  N  ,) 

1  R  01  L  xo  l,xx  o,xx  xl 

+  r(nx0owi,x0  +  wo,x0Nxei)  +  7r(N8owi,00+wo,0eNei)]=  0 

R 

The  terms  w  and  w  represent  prebuckling  rotations 
which  add  considerable  difficulty  to  the  solution  of  these 
equations.  Most  early  work  (Refs  9;  15;  18)  in  this  area 
neglected  these  terms  in  their  analyses.  With  the  computer 
available,  these  terms  do  not  have  to  be  neglected.  How¬ 
ever,  a  recent  work  by  Brush  (Ref  13)  shows  that  at  least 
for  uniform  lateral  pressure  the  effect  is  very  small. 

Figure  4  shows  that  for  most  thin  cylinders  the  effect  is 
less  than  10%. 

Therefore,  neglecting  the  prebuckling  rotations  and 
substituting  in  the  revised  kinematic  and  constitutive 
relations  we  get  a  coupled  set  of  three  linear  equations  in 
the  variables  ,  v^,  and  w^  The  solution  to  which  is 


generally  called  the  "buckling  mode." 


These  are  the  Donnell  stability  equations  in  coupled  form 
(Ref  16) .  They  can  be  partially  uncoupled  to  get  the  form 
which  are  homogeneous  in  w, .  Using  the  partial  derivative 


notation  ^  rather  than  the  indicial  notation  and  dropping 
the  subscript,  we  get: 

Vl*u  =  ZH  d3  w  +  _1_  9  3  w 
R  9x3  R3  3 0  2  3x 

_jjw_  +  _1_  9_3w  (1 

R2  303x2  R4  3 0  3 


s..  1~v  c  9  w  _  v1* 

2  3x 4 


D  v  w  + 


(N  +  |  N  Q 

X  3x2  R  x0 


3  2 


- —  +  —  =  0 

30  3x  R2  y  3  0  2 


It  is  Batdorf's  (Ref  9)  modification  of  these  equations  which 
Seide  and  Weingarten  (Ref  26)  used  to  study  an  isotropic 
cylinder  under  pure  bending  loads.  Until  the  publication  of 
(Ref  26)  the  prevailing  method  for  figuring  the  buckling 
stress  was  by  stating  that  ob  =  1.3oc  where  c?b  is  the  maxi¬ 
mum  bending  stress  and  <^c  the  critical  axial  stress  for  a 
corresponding  axial  compressive  buckling  problem.  As  Seide 
and  Weingarten  point  out,  the  value  1.3  arose  from  Flugge's 
analysis  where  he  expressed  in  terms  of  a  dimensionless 
wavelength  parameter  X  =  miTR/L;  m  denoting  the  number  of 
axial  half  waves  in  the  assumed  sinusoidal  buckle  pattern. 

He  found  ab  by  minimizing  cb  with  respect  to  X .  In  doing 
so,  he  chose  1  for  the  value  of  X  to  illustrate  the  technique, 
which  led  to  ob  =  1.3oc.  Thus  the  value  of  1.3  became  part 
of  engineering  common  practice.  Seide  and  Weingarten  used 
the  Galerkin  method  to  derive  the  stability  criterion.  They 
assumed  the  radial  deflection 
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CO 

w  =  sin  z  a  cos  n6  (13) 

L  n=0  n 

plugged  it  back  into  the  modified  Donnell  equation  and 

found  the  eigenvalue  by  matrix  iteration  of  the  determinant 

of  coefficients  for  a  .  Their  analysis  showed  that  for  all 

n 

practical  purposes  ~  oc  for  thin  cylinders.  Lakshmikan- 
tham  (Ref  18)  also  used  the  Galerkin  approach.  He  approxi¬ 
mated  the  buckle  pattern  by 

w  =  A  sin  ( x / X * )  cos  N6  for  |  0  |  =  <_  tt/2N 

=  0  elsewhere  on  the  circumference  (14) 

which  denotes  a  single  lobe  of  variable  width  dependent  on 
the  value  of  N.  He  showed  that,  with  the  single  dimple  as 
the  lowest  energy  state  z  oc  as  above.  This  method  also 
gave  an  upper  bound  =  77 when  there  exists  many 

dimples  along  the  top  half  of  the  cylinder.  They  were  also 
able  to  establish  that  generally,  the  axial  wavelength  of 
pure  bending  is  less  than  corresponding  compressive  case 
X*  £  X  with  the  equality  corresponding  to  the  axisymmetric 
case.  Thus  the  wavelength  as  well  as  the  buckle  pattern 
changes  in  going  from  the  axial  compression  to  bending 
problem.  The  wavelength  is  given  by 

X,  ~  X  =  — 2tt  ^Rt t-  for  (R/t  large)  (15) 

“  C  [  12 ( 1-v 2 ) ] 4 

=  3.45  (Rt)*5  when  v  =  .3 
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It  can  be  seen  that  this  gives  a  wavelength  quite  small  in 
comparison  with  the  length  of  the  shell. 

£  =  3.45  (|)  (|)  «  1  (R/t  large)  (16) 

More  recent  studies  have  also  used  the  Galerkin  approach. 
Reedy  (Refs  22;  23)  has  used  nondimensionalized  parameters 
in  the  Batdorf  modified  Donnell  equation.  He  uses  e  ®r' 
for  the  circumferential  modal  form  which  allows  him  to 
match,  for  sufficiently  small  values  of n  =  R9/L,  an  axial 
distribution  of  load  which  approximates  cos  6. 

The  previous  analyses  have  concentrated  on  thin, 
short  to  moderate  length  shells.  As  the  cylindrical  shell 
becomes  longer,  or  for  moderate  length  cylinders  with  ends 
that  are  not  held  circular,  nonlinear  prebuckling  effects 
must  be  considered.  Brazier  (Ref  11)  first  studied  this 
effect  in  1927  (Ref  Fig  5) .  He  noticed  that  as  an 
increasing  moment  is  applied,  the  cross-section  flattens 
and  its  moment  of  inertia  decreases  causing  ovalization. 

Due  to  the  curvature,  the  axial  stress  in  an  arbi¬ 
trary  axial  fiber  will  have  a  component  (resultant)  directed 
towards  the  neutral  axias  of  bending.  This  happens  on  both 
the  tension  and  compression  sides.  Sobel  (Ref  31)  lists 
the  following  consequences  as  a  result  of  the  oval  shape. 

1.  A  redistribution  of  stress  occurs. 

2.  The  circumferential  radius  of  curvature 


increases . 


3.  The  bi'tulinq  st  i  f  fness  of  the  cron:;-  Moot  ion  is 
reduced  as  the  moment  increases. 

4.  Nonlinear  collapse  occurs  at  a  moment  lower  than 
the  one  corespond i nq  to  bifurcation  buck  line  from  the  linear 
membrane  prebuckled  state.  Two  types  of  nonlinear  collapse 
modes  are  possible;  the  Bray,  i  or  mode  and  the  nonlinear  wavy 
mode . 

Brazier  (Kef  1 J )  determined  the  collapse  mode  by  relatinq 
the  strain  enemy  to  the  applied  moment  through  the  curva¬ 
ture  M  where  c  =  l  /,< .  His  results  qave  a  maximum 

radial  displacement . 

w  K  (17) 

under  an  appl ied  moment 


for  v  .3 


Boundary  conditions  usually  restrict  the  prebuckling  defor¬ 
mations  so  that  cross-sections  at  the  shell  edges  remain 
circular.  As  a  result,  the  flattening  occurs  at  a  higher 
stress  than  predicted  by  Brazier.  However,  as  the  L/R 
ratio  increases,  the  influence  of  the  flattening  modes 
becomes  dominant  and  the  critical  stresses  approache  that 
predicted  by  Brazier.  Figure  6  (Refs  7;  31)  details  the 
relationship  between  the  applied  moment  and  the  cross- 
sectional  flattening. 


k 

i'riiMrc:  ft.  Load  DofJection  Curves  From  ( n ) 

Akselr  d  (Ref  1)  studied  the  Brazier  effect  to  determine  a 
relevant  range  over  which  the  effect  applies.  An  asymptotic 
method  was  used  for  the  solution  of  the  prebuckling  problem 
while  an  equivalent  cylinder  approach  was  used  to  obtain 
the  critical  stress.  He  reported  that  varies  from 
.6Et/R  for  short  cylinders  to 
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—mmmm BP 


=  .295  Et/R  (R/t  large) 
~  1/2  a 

cr 

classical 


(21) 


for  cylinders  with 

L/R  2.5  (R/t)55  (22) 

His  results  are  plotted  against  discrete  points  found  by 
(Ref  7)  in  Pig  7. 


M 

c-r 

7T  Lit  2l< 


Figure  7.  Critical  Load  Factors  for 
Cylinders  Under  Bending 
'■'rom  (7) 

As  was  pointed  out  earlier,  assuming  the  same 
boundary  conditions  are  used,  test  results  generally  show 
a  lower  critical  bending  stress  than  theoretical,  mostly 
due  to  imperfections  within  the  testing  specimens.  And 
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since  the  cylinder  under  bending  has  less  of  its  surface 
at  the  peak  compression  zone,  when  compared  to  axial  com¬ 
pression,  it  is  reasonable  to  expect  that  the  maximum  bend¬ 
ing  stress  is  higher  than  the  critical  axial  stress.  Sobel 
(Ref  31)  outlines  a  method  using  a  knock-down  factor,  which 
is  an  empirical  reduction  factor,  that  accounts  for  the  dis¬ 
crepancy  between  classical  predictions  and  experimental 
results.  It  is  defined  as 

a 

C3T 

experimental 

b'c  °cr 

classical 

and  applies  to  both  bending  and  compression  in  the  same 
fashion,  respectively,  but  with  different  numerical  values. 
If  the  ratio  of  bending  to  axial  critical  stress  is  taken, 
the  following  result  is  given 


1- . 731 ( 1-e  *) 
l-.ma-e"4’) 


where  <)>  =  yg  [R/t]35 


(24) 


The  picture,  as  to  the  buckling  characteristics,  changes 
greatly  when  internal  and/or  external  stiffeners  are  added. 


Stiffening  and  Cutouts 

Stringers  and  rings  are  added  to  a  structure  to  give 
rigidity  and  reduce  imperfection  sensitivity.  If  the 
stringers  are  closely  spaced  around  the  surface,  the 
stiffened  shell  will  fail  through  general  instability.  A 
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smeared  stiffener  theory  can  be  used  to  analyze  the  total 
structure.  Its  assumptions  (Refs  21;  27)  are: 

a.  The  stiffeners  are  distributed  over  the  whole 
surface  of  the  shell. 

b.  The  normal  strains  vary  linearly  in  the 
stiffener  as  well  as  in  the  shell.  The  normal  strains  in 
the  stiffener  and  in  the  shell  are  equal  at  their  point  of 
contact. 

c.  The  shear  membrane  force  Nx^.  is  carried 
entirely  by  the  shell. 

d.  The  torsional  rigidity  of  the  stiffener  cross- 
section  is  added  to  the  shell. 

However,  if  the  stiffeners  are  spaced  relatively  far  apart 
the  use  of  discrete  theory  is  more  appropriate.  I  will  use 
the  smeared  theory  option  of  STAGSC  for  my  work.  Singer 
et  al.  (Ref  28)  did  a  parametric  study  via  a  Galerkin  sub¬ 
stitution  in  the  Donnell  stability  equations.  The  results 
showed  that  stiffener  eccentricity  has  a  large  effect  on 
axially  compressed  cylinders.  They  classified  the  effects 
as 

a.  Primary  effect-outside  stringers  increase  the 
actual  bending  stiffness  in  the  longitudinal  direction 
more  than  inside  stringers. 

b.  Secondary  effect-inside  stringers  increase  the 
actual  extensional  stiffness  in  the  circumferential  direc¬ 
tion  more  than  the  outside  stringers. 
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Figure  8  shows  the  primary  effects  whereas  Figure  9  denotes 

the  secondary  effects.  The  eccentricity  will  either  help 

or  hinder  a  restoring  moment  when  an  initial  displacement 

occurs  in  the  radial  direction.  When  the  primary  resistance 

to  buckling  is  in  the  axial  direction  the  restoring  moment, 

M  ,  remains  constant  regardless  of  the  stiffener  position. 

An  additional  moment  is  generated  by  the  stiffeners  with 

their  resultant  compressive  force  N  acting  through  the 

o 

distance  z  (the  eccentricity--also  given  as  e^  on  Fig  16) . 
Thus,  the  actual  moment  produced  by  the  shell-stiffener  com¬ 
bination  is  different  than  the  generated  shell  moment.  As 
for  the  secondary  effects,  they  come  about  due  to  an  addi¬ 
tional  membrane  force  N  produced  in  the  hoop  direction 

yo 

through  Poisson's  ratio  acting  on  the  additional  axial 

strain  caused  by  M  in  .  The  radial  component  of  this  N 

xo  yo 

will  again  assist  or  hinder  buckling  depending  on  stiffener 

location. 

The  preceding  effects  assume  the  radial  displace¬ 
ments  are  either  barrel  shaped  or  continuously  concave. 
Nelson  (Ref  21)  pointed  out  that  stringer  stiffened  pre¬ 
buckling  axial  wave  patterns  vary  from  concave  to  convex 
and  hence  should  appear  in  the  prebuckling  analysis. 

STAGSC  includes  these  effects  in  its  linear  prebuckling 
analysis. 

Another  aspect  is  where  the  load  is  applied  in  the 
shell-stiffener  combination.  Figure  10  from  (Ref  6)  depicts 
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POSITIVE  WAVE -outside  stringers 


NEGATIVE  WAVE -inside  stringers 


Figure 


NEGATIVE  WAVE -outside  stringer's 


9.  Secondary  Eccentricity  Effect 
From  (28) 
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the  results  of  loading  through  the  neutral  surface  versus 
loading  at  the  centroid  of  the  stiffeners. 


0  10  20  30  40 


Arc  Lonqth  From  Midlonqth 
in . 

Figure  10.  Prebucklina  and  Buckling 
Displacements  of  an  Internally 
Stiffened  Simply  Supported  Cyl. 

From  (6) 

/• 

It  is  apparent  from  the  Fig  10  that  the  loading 
placement  is  very  important  for  the  prebuckling  displace¬ 
ments  but  doesn't  have  too  much  of  an  effect  on  the  buckling 
mode.  Boundary  conditions  have  more  of  an  impact  there. 

Lakshmikantham  et  al.  (Ref  19)  studied  the  bending 
of  orthotropic  shells  via  the  Galerkin  procedure,  using  the 
same  assumed  deflection  function  as  for  the  isotropic 
cylinder.  They  found  ~  oc  as  before  where 
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where  Z*  =  tt 2  /  ( 1 2  ( 1-v 2 )  J*5  (26) 

anc^  Z  =  ( 1- v 2 )  ^  is  the  Batdorf  parameter.  (27) 

Referring  to  (Refs  27;  28)  further  comments  are  possible. 
Singer  noticed  that  for  very  small  values  of  Z,  that  is 
less  than  15,  inside  stiffened  shells  have  analytical 
buckling  loads  greater  than  outside  stiffened  shells.  Also, 
a  decrease  in  eccentricity  effect  occurs  for  long  shells, 

Z  greater  than  1000,  where  the  cylinders  act  unstiffened. 
Also,  contrasted  against  Z  in  the  study  were  the  stiffener 
parameters.  It  was  noted  that  as  the  stiffener  parameters 
increased,  there  was  a  corresponding  increase  in  the 
eccentric  stiffening  except  for  very  long  shells.  The 
results  showed  that  the  eccentricity  affect  behavior  depends 
very  heavily  on  the  geometry  of  the  shell;  whereas,  the 
stiffener  geometry  influences  the  magnitude  of  the  eccen¬ 
tricity  effect. 

Along  with  the  stiffener  versus  Z  studies,  the 
smeared  stiffener  theory  was  better  defined.  It  was  found 
that  the  nondimen sional  area  ratio  A^/b^t  from  Fig  16  is 
the  important  geometric  parameter  which  determines  whether 
smeared  theory  is  suitable.  Singer  (Ref  27)  concluded  that 
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with  a  stiffener  area  ratio  greater  than  0.2  stringer 
stiffened  shells  could  be  analyzed  closely  with  smeared 
theory . 

Much  study  has  also  been  directed  towards  the 
effects  of  different  boundary  conditions.  Nelson's  (Ref  21) 
thesis  was  devoted  to  a  parameter  study  using  the  standard 
delineation  of  boundary  constraints: 
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W,  =0 
'xx 
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N  =0 
xy 
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u  =  0 
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Some  of  his 

results 

were 

a.  For  all  the  simply  supported  and  clamped 
cylinders  with  and  without  cutouts,  the  u=0  axial  restraint 
is  the  most  important  factor  in  buckling  of  axially  com¬ 
pressed  cylinders.  For  small  values  of  L/R,  the  v=0  con¬ 
dition  gains  importance  in  simply  supported  stiffened 
cylinders. 

b.  For  simply  supported  cylinders,  with  and  with¬ 
out  cutouts,  internally  stiffened  cylinders  are  more  sensi¬ 
tive  to  inplane  boundary  conditions. 
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c.  The  boundary  conditions  actually  realized  in 
experimental  tests  can  be  determined  by  correlating  the 
experimental  loads  to  the  theoretical  values  for  various 
boundary  conditions.  An  agreement  between  results  would 
then  indicate  the  boundary  conditions  actually  experienced 
in  experimental  tests. 

Noting  the  latter  conclusion,  the  differences  between  experi¬ 
mental  tests  and  analytical  studies  arise  from  three  areas: 
initial  imperfections,  prebuckling  rotation,  and  boundary 
conditions.  With  judicious  forethought  the  last  two  can 
be  modeled  closely  whereas  the  first  requires  empirical 
studies . 

Theoretical  and  experimental  studies  of  cutout 
effects  under  bending  are  very  sparse,  if  any.  What  studies 
have  been  accomplished  have  dealt  mainly  with  cutout  effects 
for  axially  compressed  cylinders.  Brogan  and  Almroth  (Ref 
12)  studied  an  unstif c  ^ned  shell  with  two  symmetrical  cut¬ 
outs.  Their  theory  indicated  a  collapse  load,  30%  above  the 
experimental  load  of  the  shell  with  reinforced  cutouts, 
was  only  slightly  below  the  critical  load  for  a  cylinder 
without  cutouts.  When  they  compared  their  results  against 
Tennyson's  cylinders  with  circular  holes,  it  was  noticed 
that  the  rectangular  holes  have  less  of  an  effect.  However, 
Brush  and  Almroth  (Ref  14)  have  lately  shown  that  the  effect 
is  approximately  the  same.  Palazotto  (Ref  22)  carried  out 
a  much  more  thorough  investigation  of  cutouts  using  the 
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P^AGSA  computer  program.  He  showed  that  discrete  and 
smeared  theories  give  close  results  for  a  nodal  arrangement 
which  includes  a  circumferential  mesh  line  at  and  between 
stringers  and  axial  mesh  spacing  equal  to  the  stringer 
spacing.  This  arrangement  produced  results  within  3%  of 
those  with  a  more  detailed  mesh  arrangement.  Further,  a 
mesh  arrangement  that  insures  at  least  five  nodes  per  half 
sine  wave  was  shown  to  be  sufficient.  Nelson  (Ref  21) 
further  studied  the  effects  of  sn  eared  stringer  reinforce¬ 
ment  of  the  cutouts.  He  found  that  cutouts  affected  the 
externally  stiffened  shells  to  a  greater  extent  than 
internally  stiffened  shells  and  this  effect  diminished  as 
the  L/R  ratio  increased.  Also,  he  found  that  the  influence 
of  cutouts  on  the  critical  loads  of  stringer  stiffened 
cylinders  is  analogous  to  the  effects  that  initial  geometric 
imperfections  have  on  theoretical  loads  for  "perfect" 
cylinders. 

This  introduction,  although  not  complete  as  to  all 
the  studies  which  have  been  completed  in  this  area,  should 
serve  as  an  overview  of  the  subject. 

General  Procedure 

The  Skogh-Brogan  bending  moment  routine  (Ref  29) 
is  combined  with  the  STAGSC  (Ref  5)  computer  code  to  analyze 
the  effects  on  cylindrical  cylinders.  The  L/R  ratio  is  kept 
very  short  to  minimize  the  Brazier  effect.  Using  the  linear 
bifurcation  branch,  a  rotation  is  applied  to  cylinders 
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without  stiffening.  The  impetus  is  to  study  the  deflection 
buckling  data.  The  cylinder  L/R  is  then  increased  up  to  20 
to  1  to  see  how  the  bending  boundary  deflections  are  changed. 
Then,  end  rings  are  added  to  the  shell  to  see  how  the  near 
bending  boundary  deflections  are  effected.  The  last  half 
of  the  study  has  strinqers  added  to  the  cylinder  for 
analysis,  with  and  without  12"  cutouts.  The  STAGSC  smeared 
stiffener  theory  is  used  for  this  analysis.  A  parameter 
study  is  accomplished  with  the  L/R  ratio  changing  from 
.698  to  1.361. 
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II.  STAGS  Procedures 


STAC', SC  Theory 

References;  2,  3,  4,  5,  6,  and  13  contain  to  varyinq 
degrees  the  basis  on  which  the  STAGSC  code  is  bui  it.  I 
recommend  (Refs  3;  4)  as  a  primer  before  employing  the 
users  manual  (Ref  5),  since  the  proposed  theory  section  for 
the  STAGSC  manual  was  never  published.  Also,  a  newer  ver¬ 
sion  STAGSC- 1  is  now  out  in  the  field  which  mainly  uses  a 
finite  clement  approach.  Therefore,  it  will  probably  be 
harder  in  the  future  to  got  help  with  the  STAGSC  version. 

One  advantage  that  the  C-l  version  should  have,  although  it 
is  not  available  yet,  is  a  post-processing  graphics  package. 
For  any  structural  code'  a  picture  representation  is  very 
much  needed  to  get  a  grasp  on  what  the  surface  is  doing. 
Without  it,  the  time  devoted  to  plotting  points  makes  that 
particular  code  less  worthwhile.  I  will  outline  the  basics 
of  the  STAGSC  finite  difference  process  in  finding  the 
buckling  values,  then  indicate  the  major  differences  between 
STAGS  A  and  C  versions  for  finite  difference  analysis,  and 
mention  some  of  the  facets  involved  in  its  smeared  theory. 

Structural  behavior  is  governed  by  the  equation 

Mil  +  D(u)  n(u)  +  K  (u)  F  (27) 
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where  u  is  a  vector  of  displacement  components,  M  is  the 
mass  matrix,  D  and  B  include  forces  that  are  functions  of 
deformation  velocity  and  structural  deformation  respec¬ 
tively,  K  is  the  generally  nonlinear  stiffness  operator,  and 
F  a  vector  of  external  forces.  If  the  loads  vary  slowly 
all  time  derivatives  can  be  neglected.  Figure  11  comes  from 
(Ref  4) .  It  delineates  the  solution  avenues  available  for 
the  different  applicable  terms  of  Eq  27.  Assuming  the 
external  forces  vary  very  slowly  and  the  load  is  independent 
of  deformation,  then  STAGSC  uses  a  potential  energy  method 
like  Eq  6  to  develop  the  applicable  equations.  A  surface 
mesh  is  introduced  and  derivatives  of  the  displacement  com¬ 
ponents  are  expressed,  using  a  truncated  Taylor  series 
expansion,  in  terms  of  the  displacements  at  the  mesh  points 
by  different  approximations.  Brush  (Ref  14)  gives  a  good 
introduction  to  the  subject  and  solves  an  example  column 
buckling  problem.  The  derivative  approximations  are  substi¬ 
tuted  into  the  strain  energy  equations.  For  stability,  it 
was  related  earlier  that  the  second  variation  of  V  is  set 
equal  to  zero.  Vanishing  of  the  first  variation  of  V  leads 
to  the  equation 

LX  =  F  (28) 

where  L  is  the  stiffness  operator  that  relates  the  finite 
nodal  displacement  vector  X  to  the  external  forces.  For 
the  bifurcation  problem  if  another  displacement  vector  Y  in 
the  neighborhood  of  X  satisfies  the  equation,  then 
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LY  =  F  (29) 

And  since  there  can  be  multiple  solutions  to  Eq  28,  a  neces¬ 
sary  condition  for  bifurcation  is  that  L'x  be  a  singular 
matrix  where  l!  denotes  the  derivative  of  the  stiffness 
operator.  This  leads  to 

det  (I/x)  =  0  (30) 

It  is  then  assumed  that  X  may  be  written  as 

X  =  XXL  (31) 

where  X  is  a  linear  solution  for  the  vector  FT ;  this  leads 

Li  h 

to  Eq  30  becoming 


det  (L'Axl)=0  (32) 

This  is  equivalent  to  an  algebraic  eigenvalue  problem.  Its 

solution  for  A  gives  the  critical  load  as  a  multiple  of  the 
linear  force  vector  FT 

Li 

fcr  -  Xfl  1331 

In  this  discussion,  it  was  mentioned  that  a  certain  differ¬ 

ence  scheme  is  used  to  evaluate  the  derivatives.  STAGSA 
was  based  on  a  half  station  scheme.  See  Fig  12  from  (Ref 
3) .  This  scheme  has  a  weakness  in  its  inefficiency  for 
nonlinear  or  stability  analysis.  It  also  has  the  disadvan¬ 
tage  that  in-plane  displacements  must  be  defined  at  ficti¬ 
tious  points  outside  the  shell  boundaries.  The  STAGSC 
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Fiqure  12.  STAGSA  and  STAGSC  Finite  Difference  Schemes 
From  (3) 


version  uses  the  whole  station  scheme.  The  point  in  ques¬ 
tion  is  surrounded  by  the  total  integration  area  broken 
down  into  4  smaller  regions.  The  corner  of  each  smaller 
area  has  defined  the  requisite  quantities,  including  the 
rotations.  This  method  allows  better  convergence.  It 
has  been  postulated  in  (Ref  3}  that  the  finite  difference 
approach  converges  toward  the  correct  solution  if: 

1.  The  local  truncation  error  vanishes  with  the 
gridsize. 

2.  For  small  values  of  nodal  separation  distance, 
the  solution  varies  continuously  with  the  input  data  loads. 

It  was  brought  out  in  the  background  section  that 
Singer  (Ref  27)  gave  several  conditions  for  using  the 
smeared  stiffener  analysis.  Since,  for  my  work  the  stiff¬ 
ener  area  ratio  is  .9333  >  .2  recommended  by  Singer,  the 
smeared  theory  is  applicable.  Brush  (Ref  14) , 
Lakshmikantham  (Ref  19),  and  Singer  (Ref  28)  outline  how 
smeared  theory  is  incorporated  into  the  equilibrium  equa¬ 
tion.  It  is  handled  through  the  generalized  constitutive 
equations  where  the  coupling  effects  are  added  to  or  sub¬ 
tracted  from  the  existing  stiffness  parameters.  As  an 
example,  the  primary  axial  compression  effect  of  a  smeared 
stiffener  is  added  to  the  membrane  extensional  stiffness 
parameter  C  expression.  It  has  the  form 


C  =  C  + 
new  old 


where  Ag,  ds  are  the  area  of  and  distance  between  stringers. 
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The  STAGSC  (Ref  5)  code  lets  the  user  detail  the 
shell  wall  properties  and  dimensions,  either  directly 
through  regular  data  cards  or  by  a  user  written  subroutine 
Wall.  The  program  further  handles  the  addition  of  smeared 
or  discrete  stiffeners  along  the  surface,  and  the  mesh  grid 
size  can  be  varied  through  specification  in  the  discretiza¬ 
tion  section.  Since  my  analysis  includes  nonsymmetric 
cutouts,  the  regular  data  cards  were  unable  to  handle  the 
specifications,  and  I  had  to  use  the  user  defined  subroutine 
Wall  to  input  my  data.  Appendix  B  contains  two  listings  of 
the  subroutine  Wall.  The  first  listing  does  not  include 
cutout  data  whereas  the  second  listing  does.  Although  for 
the  case  of  no  cutouts  (listing  1)  the  use  of  the  Wall  sub¬ 
routine  was  not  needed,  I  decided  to  use  it  for  uniformity 
so  that  the  cutout  results  could  be  compared  directly.  In 
both  cases  smeared  theory  was  incorporated  within  the  Wall 
subroutine.  The  mesh  dimensions,  surface  parameters,  stiff¬ 
ener  configuration,  and  cutout  relationships  will  be  out¬ 
lined  in  the  modeling  section  of  this  report. 

The  Skogh-Brogan  Routine 

To  apply  pure  bending  moments,  in  previous  work  done 
with  STAGS  (Refs  5;  15;  32),  an  end  line  load  was  varied 
around  the  edge  uniformly  with  the  cosine.  (See  Fig  13.) 

N  =  -N  cos  0  (35) 


Possible  Ring  Locations 


Figure  13.  Applied  Bending  Moment  Through  Uniformlv 

Varying  Axial  Loading 


where  N  is  the  amplitude  of  the  load  and  0  is  the  circum¬ 
ferential  coordinate.  The  generated  bending  moment  is 
therefore 


M  =  NirR2  (36) 

This  method  did  impart  the  desired  moment  numerical  values 
however,  it  would  also  lead  to  some  warping  of  the  end  plane 
resulting  in  a  load  that  was  not  pure  planer  bending,  and 
in  some  cases  a  net  axial  force  also  developed.  End  rings 
were  used  to  try  and  counteract  the  warping  tendency. 

Mr.  Jorgen  Skogh  and  Frank  Brogan  (Ref  29)  working 
at  Lockheed's  Palo  Alto  Research  Laboratory  approached  the 
problem  from  a  different  viewpoint.  They  set  out  to  enforce 
a  planer  end  rotation  and  accept  the  applied  moment  value 
which  resulted.  They  were  able  to  accomplish  this  through 
the  use  of  constraint  equations  which  required  that  any 
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point  on  the  end  plane  meet  certain  boundary  kinematic 
restraint  relationships.  Thus  free  end  rotation  became 
the  independent  variable  (see  Fig  14) .  The  free  end  of 
the  cylinder  is  shown  displaced  through  an  angle  <P  which  is 
assumed  to  be  small  so  that  cos  4>  ~  1.  The  equations  which 
they  used  to  keep  the  end  plane  circular  are 

v  sin  9  -  w  cos  6  =0  (37) 

-v  +  v  cos  8  +  w  sin  0=0  .  (38) 

o 

Eqs  37  and  38  keep  the  triangle  formed  by  v,  w,  and  vq 
(Fig  14)  in  equilibrium.  That  is  why  the  angle  4>  should  be 
kept  small — to  keep  the  length  vq  linear  in  Eq  38.  The 
requirement  that  this  circular  end  remain  planer  under  rota¬ 
tion  is  fulfilled  by 

uo  =  (Ugo-U0)  sin  0  -  U  =  0  (39) 

This  says  that  for  any  angle  9  (position  around  the  circle) 
the  vertical  distance  u  is  a  function  of  the  difference 
u90  “  u0  as  in  14.  By  changing  that  difference, 

the  distance  u  varies,  changes  for  a  given  uQ  and  0.  In 
effect,  Ugg  -  uQ  is  a  measure  of  the  slope  which  has  changed, 
meaning  the  angle  of  rotation  has  changed. 

Skogh  and  Brogan's  next  ingenious  step  was  to  input 
this  difference  as  a  displacement  on  a  "dummy"  branch.  The 
term  dummy  is  used  in  the  respect  that  this  second  branch 
is  only  used  as  a  method  for  inputing  data  and  has  no  other 
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relevance.  STAGSC  has  the  ability  to  match  boundaries 
along  many  connecting  branches.  In  this  case  the  boundary 
is  not  connected.  However,  the  second' branch  is  used  to 
input  the  desired  displacement.  Skogh  reads  in  the  desired 
displacement  u^  via  load  data  cards  for  the  second  branch. 
He  then  substitutes 


u 


90 


-  u  =  u. 


(40) 


into  Eq  30;  and  along  with  the  other  two  constraint  Eqs  37, 
38  forms  the  subroutine  UCONST.  UCONST  allows  the  STAGS 
user  to  provide  kinematic  constraints  in  the  displacement 
unknowns  at  internal  and  boundary  points.  The  general  con¬ 
straint  relation  is  of  the  form 


N 

I  C.V.  =  0  (41) 

i  11 

where  are  constants  and  are  the  displacement  com¬ 
ponents.  What  happens  internally  is  that  when  the  equa¬ 
tions  are  being  set  up  for  the  individual  mesh  points,  at 
those  points  where  UCONST  is  applicable  Eq  41  has  to  be 
solved  for  also.  The  displacements  are  normally  not  fixed 
but  rather  have  to  meet  certain  kinematic  relationships  with 
surrounding  points  while  the  entire  structure  satisfies  com¬ 
patibility  and  equilibrium.  This  is  the  difference  between 
using  UCONST  vs  input  boundary  constraint  cards,  which  form 
fixed  positions.  The  dummy  branch  is  used  to  internally 
impart  the  desired  rotation.  Therefore,  by  specifying  the 
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top  row  as  the  place  of  interest  in  UCONST,  the  desired  con¬ 
straint  equations  are  inforced  with  -  uq  =  u^.  The 

UCONST  subroutine  is  listed  in  Appendix  B. 

A  few  points  about  the  actual  operation  of  the 

method. 

1.  The  second  dummy  branch  can  be  of  any  type,  so 
the  easier  to  model  the  better.  I  chose  a  1  x  1  inch  plate. 

2.  The  specified  displacement  on  this  branch  is  in 
the  x  direction.  Changing  the  sign  on  the  input  reverses 
the  slope  in  UCONST  and  hence  the  applied  rotation. 

3.  When  trying  to  run  a  bifurcation  analysis  of 
the  primary  branch,  STAGS  will  automatically  run  a  bifurca¬ 
tion  analysis  of  all  branches  and  stop  when  the  first  con¬ 
vergence  criteria  is  met  for  any  branch;  not  all  branches. 
Hence,  to  make  sure  that  the  bifurcation  analysis  is  car¬ 
ried  out  on  the  branch  of  interest,  the  other  branches 
should  have  their  bifurcation  analysis  suppressed.  This  is 
accomplished  by  setting  the  variable  L1N=1  on  the  II  card 
for  all  the  nondesired  branches.  See  Appendix  B1  for  an 
input  card  listing  and  the  user's  manual  (Ref  5)  for  further 
guidance. 

4.  UCONST  requires  a  do  loop  to  set  its  way  through 
all  the  mesh  points.  When  the  nodal  size  is  changed  in  the 
circumferential  mesh  direction:  the  dimension,  data,  do 
loop,  singularity  check,  and  IY  cards  have  to  be  updated 

in  the  UCONST  subroutine  as  well  as  the  number  of  constraint 
relations  variable  NCONST  on  the  B1  card. 


III.  Model in 


Shell  Model 

Figure  15  denotes  the  pure  bending  rotation  conven 
tions.  The  -90°  fibres  are  in  a  tension  field  whereas  the 
+90°  fibres  are  in  a  compression  field.  The  radius  R  is 
57.2958  inches.  This  facilitates  measurements,  since  one 
inch  on  the  circumference  equates  to  one  degree.  The  cut¬ 
outs  are  dashed  indicating  that  they  are  not  always  there 
in  my  analyses.  When  they  are,  2a  =  12  inches. 


-90° 

Figure  15.  Pure  Bending  Rotation 


Figure  16  is  another  view  of  the  cylindrical  shell 
In  addition,  it  has  the  stringer  data  drawn  and  tabulated. 
Palazotto  (Ref  22)  demonstrated  nonlinearities  due  to  cut¬ 
outs  in  stiffened  shells  were  small  even  for  a  cutout  size 
of  2a  equal  to  .3L.  To  keep  w  chin  that  finding  Nelson 
(Ref  21)  used  smaller  cutouts  when  he  studied  shorter 


R=57 . 2958 " 
t=.  1" 
2a=12" 
b,=3" 
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Figure  16.  Shell,  Strinner  Geometry  and  Sion  Conv 


cylinders.  Since  the  major  part  of  my  work  will  be  with  an 
L/R  '  1  range,  I  will  also  use  the  smaller  cutouts. 

The  stringers  placement  followed  that  of  Palazotto 
and  Nelson.  The  three  inch  on  center  spacing  with  Aj=.28 
is  well  within  the  realm  of  smeared  stiffener  theory. 

A  significant  effect  of  studying  bending  problems 
is  that  the  moment  takes  away  the  third  plane  of  symmetry 
due  to  loading  asymmetry.  It  is  not  axi symmetric .  So  a 
quarter  of  the  shell  has  to  be  analyzed  each  run.  This 
doubles  the  nodal  requirements  for  the  same  "fineness"  of 
spacing  which  modeling  one  eighth  of  the  cylinder  allows 
for  axial  compression  analysis.  However,  I  am  therefore 
able  to  study  asymmetrical  cutouts  (one  rather  than  two) 
without  changing  the  modeling  strategy.  The  boundary  con¬ 
ditions  are  defined  in  Fig  17.  Conditions  of  symmetry  are 
imposed  on  boundaries  2,  3  and  4  throughout  the  analysis. 

The  simple  support  boundary  conditions,  Kqs  26,  presented 
earlier  are  imposed  on  the  shell  ends  acting  under  axial  com¬ 
pression  along  with  the  free  condition  for  bending  analysis. 
For  axial  compression  runs,  both  line  loads  and  displace¬ 
ments  were  imposed  on  the  shell  ends.  For  bending  analysis, 
the  dummy  branch  induced  rotation  is  always  used. 

When  building  some  runs  for  comparison  with  other 
published  results,  wo  noticed  what  has  turned  out  to  be  an 
internal  computational  error  in  the  STAGSC  program.  Figure 
18  shows  STAGS  results  for  axial  compression  runs  under  dis¬ 
placement  and  force  loading  contrasted  against  Nelson's 
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BOUNDARY 
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PREBUCKLING 
Free,  SSI 
Symmetry 
Symmetry 
Symmetry 


BUCKLING 
Free,  SS3 ,  SS4 
Symmetry 
Svmmetry 
Symmetry 


Note:  Symmetry  means  u  =  Rj=  0  on  lines  1  and  3 

v  =  0  on  lines  2  and  4 


Figure  17.  Boundary  Conditions  for  the  Shell  Model 
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Figure  18.  Comparison  of  Axial  Compressive 
Failure  Loads  Under  Force  (SS4) 
and  Displacement  (SS3)  Conditions 


results.  Data  for  stringer  stiffened  cylindrical  shells 
with  and  without  cutouts  are  examined.  All  Nelson's  runs 
were  made  with  displacement  loading  while  the  boundary 
conditions  varied.  The  comparative  results  for  SS4  condi¬ 
tions  are  reasonably  close.  It  is  when  SS3  conditions  are 
compared  that  a  large  disparage  exists:  around  450,000  for 
Nelson's  vs  800,000  for  present  work.  However,  when  force 
loading  is  introduced,  A  value  for  the  SS3  cutout  run,  the 
results  fall  within  1%  of  each  other.  This  has  led  to  the 
conclusion  that  when  the  displacement  loading  is  used  with 
the  STAGSC  bifurcation  option  the  stiffness  matrix  formed 
for  the  prebuckling  boundary  conditions  is  not  decoupled 
from  the  buckling  boundary  conditions.  When  force  loading 
is  used  this  does  not  happen  since  a  specific  u  has  not  been 
specified.  When  the  force  loading  runs  are  made  with  SS4 
boundary  conditions,  the  same  results  are  obtained  as  with 
the  SS3  displacement  combination.  The  first  clue  to  the 
reason  for  these  results  came  from  the  eigenvector  printout 
sheet.  For  the  displacement  runs  with  SS3  conditions, 
the  top  row  u  values  were  zero  which  they  shouldn't  have 
been.  This  did  not  show  up  when  force  loading  was  used. 

Mr.  Jorgen  Skogh,  Lockheed's  STAGSC  consultant  and  monitor, 
has  been  contacted  about  these  findings. 

This  development  has  caused  us  to  compare  the  bend¬ 
ing  moment  results  against  the  force  SS4  combination  axial 
compression  results.  Although  the  top  row,  boundary  one 
condition,  is  set  as  free  for  the  bending  runs,  we  feel  that 
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by  setting  in  a  definite  rotation  angle,  the  same  thing 
is  happening  with  the  stiffness  matrix.  Results  have  shown 
this  to  be  the  case.  Future  work  in  this  area  with  STAGSC 
should  center  around  resolving  this  internal  mathematical 
programming  deficiency. 

Finite  Difference  Mesh  Arrangements 

The  grid  systems  developed  for  shell  models  with  and 
without  cutouts  are  illustrated  in  Appendix  C.  Since 
smeared  theory  is  used  in  this  study,  a  mesh  line  should  be 
placed  under  each  stringer  in  the  circumferential  direction. 
This  necessitates  a  minimum  of  61  nodes  in  the  circumferen¬ 
tial  direction.  Convergence  studies  were  made  with  the 
results  confirming  the  need  for  this  arrangement  in  the 
compressive  zone.  However,  there  seems  to  be  room  for 
economy  in  the  tension  zone.  A  comparative  bending  run 
(Fig  19)  with  half  the  number  of  mesh  lines,  that  is  6  inches 
apart,  on  the  tension  side  was  made  with  the  same  spacing  as 
before  on  the  compressive  side.  The  buckling  loads  computed 
out  to  be  the  same  but  the  eigenvectors  were  slightly  dif¬ 
ferent.  More  importantly,  this  was  done  on  an  unstiffened 
cylinder.  Future  studies  will  have  to  determine  whether 
economies  can  be  gained  using  this  method.  For  this  study, 
the  mesh  arrangement  has  remained  symmetric  about  the  middle 
column.  Another  point  should  be  made  here.  Previous 
studies  using  STAGSA  (Refs  15,  21,  22,  32)  have  used  as  a 
rule  of  thumb,  5  nodes  per  half  sine  wave  in  the  eigen 
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analysis  as  a  good  measure  for  convergence.  As  was  brought 
out  earlier,  STAGSA  used  a  half  station  difference  scheme 
whereas  STAGSC  uses  a  whole  station  method.  Almroth 
(Ref  3)  did  converge  comparisons  between  the  two  programs 
showing  STAGSC  to  be  much  quicker  and  accurate.  In  fact, 
for  similar  accuracy  we  can  use  almost  twice  as  large 
(sparser)  grid  spacing.  This  indicates  that  the  smaller 
node  count  per  half  sine  wave  is  allowable.  However,  you 
will  want  to  be  able  to  model  the  sinusoidal  behavior 
accurately. 

Since  this  analysis  is  concerned  with  the  displace¬ 
ment  effects  of  applying  a  pure  rotation  to  the  end  of 
cylinder,  an  axial  mesh  spacing  of  1"  near  the  top  boundary 
is  used  for  all  runs.  The  larger  pertinent  shape  changes 
are  seen  to  take  place  within  the  top  5  or  6  inches  as  was 
also  indicated  by  Nelson  (Ref  21) .  Stiffened  cylinders 
generally  buckle  with  one  half  sine  wave  in  the  axial  direc¬ 
tion.  So  spacing  is  determined  more  by  economy  and  the 
desire  to  keep  the  mesh  as  square  as  possible.  The  3" 
axial  mesh  spacing  beneath  the  top  5  or  6  rows  is  entirely 
adequate.  Unstiffened  cylinders  display  more  waves  in  the 
axial  direction  as  discussed  under  background.  The  3" 
axial  spacing  gives  2h  nodes  per  half  sine  wave  axially 
which  I  believe  is  sufficient  for  STAGSC  analysis  if  the 
circumferential  mesh  spacing  is  "fine"  (see  Fig  19) . 

For  the  cutout  region  Palazotto  (Ref  22)  and  Nelson 
(Ref  21)  have  shown  that  a  mesh  spacing  of  1.5"  extending 
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15"  from  the  center  of  the  cutout  represent  an  optimum 
refinement  as  far  as  computer  time  vs  accuracy  is  concerned. 

I  will  use  this  spacing  strategy,  as  it  also  gives  an  addi¬ 
tional  mesh  line  between  stringers  circumferentially. 

For  the  compressive  analyses,  I  have  used  a  com¬ 
bination  of  1/8  shell  and  1/4  shell  models.  The  axial 
dimensions  remain  the  same  as  discussed  earlier.  The  cir¬ 
cumferential  dimensions  mainly  vary  in  the  total  length  of 
the  shell  modeled;  not  the  spacing  between  the  node  points. 
That  is,  for  unstiffened  axial  compression,  I  ran  runs  with 
15  by  31  spacing  over  90°  circumferentially  about  15  by  61 
spacing  over  180°  circumferentially .  The  total  number  of 
node  points  doubled  but  the  spacing  remained  the  same  between 
them.  This  was  done  in  a  study  of  boundary  condition  effects. 
For  the  stiffened  shells  without  cutouts,  I  kept  the  circum¬ 
ferential  spacing  the  same  as  for  the  bending  analysis:  61 
columns  over  180°.  For  the  cutout  analysis,  I  used  1/8 
shell  modeling  as  did  Nelson  (Ref  21) . 
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IV.  Results 


Unstiffened  Analyses 

Short  shells  were  originally  chosen  for  our  model 
analysis  in  order  to  stay  away  from  the  Brazier  non  linear 
effects.  Also,  previous  work  had  been  done  (21,  22)  on 
stiffened  short  cylinders  with  cutouts  under  axial  com¬ 
pression  against  which  comparisons  could  be  made.  Un¬ 
stiffened  cylinders  were  first  examined  to  give  us  some  ex¬ 
perience  with  STAGSC . 

Figure  19  is  a  plot  of  stress  intensities  for  the  max¬ 
imum  compressive  zone  at  buckling  for  different  mesh  arrange 
raents .  The  reasons  for  these  grid  dimensions  were  given  in 
the  modeling  section.  Models  numbered  2,  3,  4,  6  and  7, 
shown  in  Appendix  C,  were  used  for  these  convergence  tests. 
For  axial  compression  analysis,  the  15  x  31  mesh  over  180° 
is  too  sparse.  The  buckling  values  obtained  with  it  were 
36%  of  classical.  Keeping  the  axial  mesh  the  same  and 
halving  the  circumferential  spacing  gives  results  within  78% 
of  classical.  Again  halving  the  circumferential  spacing, 
that  is  doubling  the  node  points  in  that  direction  to  121, 
the  answer  comes  to  within  95%  of  classical  values.  However 
the  total  node  point  count  of  1815  is  very  costly  to  solve. 
It  was  decided  that  the  characteristics  we  were  looking  for 
in  the  unstiffened  cylinders  could  be  modeled  adequately  by 
the  intermediate  15  by  61  mesh  arrangements.  Also,  as  is 
seen  later  on,  for  stiffened  shells  (Fig  28)  this  same  grid 
spacing  gives  very  good  results. 


Figure  19  shows  the  circumferential  and  axial  node 
point  counts  per  half  sine  wave  for  the  respective  mesh 
arrangements.  As  the  circumferential  spacing  became  "finer" 
(more  mesh  lines  per  direction)  the  node  point  count  rose 
accordingly.  For  121  mesh  lines  we  get  6  node  points  per  % 
sine  wave.  It  is  interesting  that  for  a  given  axial  spacing 
the  convergence  is  shown  to  be  a  function  of  the  circumfer¬ 
ential  spacing.  This  means  that  the  axial  spacing  has  to  be 
at  least  a  certain  minimum  to  model  the  cylinder  adequately. 
There  has  to  be  enough  grid  points  to  let  the  equations 
model  the  shells  stiffness  correctly.  However,  it  seems 
that  a  lack  of  mesh  lines  in  one  direction  can  be  compen¬ 
sated  for,  up  to  a  point,  by  increasing  the  lines  in  the 
other  direction.  The  best  practice  is  to  keep  the  grid 
sizing  as  square  as  possible  with  as  fine  a  mesh  spacing 
as  is  economically  feasible. 

Also,  mesh  spacing  strategy  usually  starts  with  a  large 
grid  spacing,  then  decreases  the  spacing  around  regions  of 
interest.  This  method  was  used  in  this  thesis,  by  having  1 
grid  line  per  inch  near  the  top  boundary  where  the  loading 
is  applied,  and  1^  inch  spacing  around  the  cutouts.  However, 
this  was  done  to  study  the  deflections  and  stresses  in  these 
regions,  not  to  get  a  better  buckling  convergence.  As  was 
mentioned  earlier,  for  bending  analyses,  only  the  side  of 
the  cylinder  under  the  compressive  zone  is  of  buckling  in¬ 
terest.  Figure  19  shows  that  a  mesh  arrangement  (model  6) 
with  half  the  mesh  lines  in  the  tension  zone,  can  give  the 
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same  results  for  bending  analyses  as  symmetrical  arrange¬ 
ments.  This  was  only  tested  for  unstiffened  cylir  !ers; 
further  studies  will  have  to  investigate  this  mesh  strategy 
for  stiffened  cylinders  where  common  usage  says  there  should 
be  at  least  one  grid  line  under  each  stiffener.  For  all 
the  axial  compression  vs  end  rotation  loading  analyses  in 
this  thesis,  the  circumferential  grid  spacing  was  symmetri¬ 
cal  about  the  middle  column. 

The  major  finding  of  this  thesis  can  be  seen  in  Figs 
20,  25,  26  and  27.  Looking  first  at  unstiffened  short  cyl¬ 
inders  (Fig  20)  under  pure  bending,  via  kinematic  end  ro¬ 
tation  constraints,  the  prebuckling  deflection  shapes  do  not 
follow  what  is  normally  called  "beam  bending"  deflection. 
That  is,  instead  of  the  cylinder  deflecting  as  pure  body 
motion  with  the  corresponding  convex  shape,  the  cylinder 
deforms  as  shown  in  Fig  20.  This  figure  details  only  the 
top  half  of  the  cylinder  with  the  bottom  half  being  a  mirror 
image.  There  isn't  enough  material  in  these  short  unstiff¬ 
ened  cylinders  to  allow  the  stress  loading  caused  by  the 
applied  end  rotation,  to  dissipate  through  flexure  type  de¬ 
flections.  Another  way  of  saying  it  is  that  the  cylinder 
isn't  stiff  enough  to  counteract  the  local  end  moments  Mx 
which  cause  a  shear  strain  type  deformation.  By  keeping 
the  end  plane  undistorted  the  material  around  -90°  is  put 
under  pure  tension;  whereas,  the  fibers  around  +90°  are 
subjected  to  pure  compression.  The  resulting  deformation 
has  the  material  on  the  tension  side  necking  in,  and  the 
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material  on  the  compression  side  bulging  out.  The  maximum 
extent  of  the  inward  deformation  or  outward  bulge  is  labeled 
MWD  (on  Fig  20) .  Also  on  this  figure,  the  V  deflection  at 
the  middle  column  is  sketched.  This  line  is  seen  to  take 
the  shape  normally  associated  with  bending.  This  indicates 
that  the  vertical  plane  through  0°  and  180°  about  which  the 
end  rotation  is  applied  does  deflect  in  the  "normal  manner". 
Figure  21  is  a  top  view  of  the  first  few  rows  of  the  cylin¬ 
der.  It  shows  how  the  deformations  vary  around  the  circum¬ 
ference  increasing  to  a  maximum  at  row  5  (4"  down).  At  0° 
and  180°  there  is  a  very  slight  motion  (not  discernible  in 
this  plot)  towards  the  bottom  of  the  page. 

It  should  be  mentioned  at  this  point  that  the  actual 
deflections  are  of  order  10~3  inch  with  an  L/R  ratio  of 
1.047.  So,  this  effect  would  not  be  discernible  to  the 
naked  eye. 

Again  looking  at  Fig  20,  the  inward/outward  deflections 
in  the  tension  and  compressive  zones  respectively  do  not  get 
back  to  the  original  circumference  position.  This  changes 
when  the  L/R  is  doubled  to  2.09  where  the  form  takes  on  the 
shape  designated  in  Figs  22  and  23.  The  length  of  cylinder 
is  doubled  keeping  the  radius  the  same  while  the  applied  end 
rotation  angle  is  kept  constant.  The  result  is  that  the  de¬ 
formations  now  cross  over  the  original  circumferential 
position  and  start  to  assume  flexure  deflection  type  displace¬ 
ments  at  the  mid  section.  This  transfer  region  is  called 
the  crossover  point  (CP)  and  is  designated  in  Fig  22.  As 
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the  length  to  radius  ratio  increases  for  the  same  applied 
input  rotation,  the  deformation  effects  diminished  in  magni¬ 
tude  as  well  as  position  along  the  length  of  the  cylinder. 
Figure  22  shows  that  as  the  L/R  ratio  approaches  20  the 
effects  disappear.  In  fact,  for  L/R  ratios  greater  than 
6  the  effects  can  be  said  to  be  negligible.  This  is  still 
well  below  the  L/R  value  of  50  where  the  Brazier  effect 
starts  to  take  place.  Figure  23  is  a  plot  of  the  W  deflec¬ 
tion  vs  V  deflection  at  the  mid  section  as  the  L/R  varies. 
Since  the  V  deflection  of  the  middle  column  has  been  shown 
to  follow  "beam  type"  flexure,  the  relation  of  the  W  de¬ 
flection  at  the  end  columns  to  V  gives  a  good  indication 
of  the  cylinders  behavior.  As  is  seen  in  Fig  23,  the  ratio 
is  almost  unity  at  L/R  equal  to  20  and  for  all  practical 
purposes  can  be  considered  as  such  for  L/R  ratios  greater 
than  6. 

A  final  comment  about  using  the  pure  bending  loading 
method;  all  motion  degrees  of  freedom  are  accounted  for 
with  respect  to  free  body  motion  through  the  symmetric 
conditions  on  boundaries  2,  3,  and  4  except  for  translation 
in  the  -90°/90°  plane.  Some  restriction  must  be  added  in 
the  program  to  stop  this  from  occuring.  A  restriction  of 
setting  V  equal  to  0  at  the  middle  column  of  row  1  was  used 
This  was  chosen,  rather  than  restricting  V  at  the  center, 
since  the  end  plane  rotation  is  the  applied  loading  method. 
The  V  restriction  at  the  center  would  lead  to  the  same  re¬ 
sults  but  the  ends  would  be  deflected  rather  than  the 
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center.  Also,  the  V  restriction  is  preferable  over  a  W  re¬ 
striction  at  the  -90°  or  90°  columns  since  V  restricts  in¬ 
plane  motion;  whereas,  the  W  would  restrict  motion  normal 
to  the  surface  which  is  much  weaker. 
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Figure  19.  Buckling  Stress  Intensities  vs  Mesh  Size 
For  Pure  Bending  and  Axial  Compression  of 
an  Unstiffened  Cylinder 
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Figure  20.  Prebuckling  Deflections  For  Unstiffened 

Cylinder 
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Compression  Zone 


Figure  21.  Top  View  of  an  Unstiffened 
Cylinder  Under  Pure  Bending 


Figure  22.  Length  to  Radius  Effects  For  an  Unstiffened 
Cylinder  With  an  Applied  End  Rotation 
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W  Deflection  at  Cylinder  Center 
Plane  ie  Bottom  row  of  Model, 
column  1 

V  Deflection  also  at  Cylinder 
Center  Plane  but  at  Middle 
Column  of  the  Model 


Since  the  unstiffened  cylinders  were  not  able  to  re¬ 
sist  the  deformations  caused  by  the  end  plane  bending  ro¬ 
tation,  a  study  using  discrete  rings  positioned  near  the 
ends  was  undertaken.  Figure  24  is  a  plot  of  the  results 
for  the  first  7  rows  of  the  compression  column  at  90°.  It 
is  seen  immediately  that  for  the  ring  configurations  cho¬ 
sen,  the  basic  deformation  pattern  does  not  change. 

One  ring  added  at  the  top  row,  whether  internal  or  ex¬ 
ternal,  has  no  effect.  It  acts  like  a  point  mass  and  does 
not  affect  the  axial  bending  stresses.  Rings  added  at  row  2 
(1" )  or  row  3  (2")  do  affect  the  initial  slope  and  specific 
shape  of  the  deformation  pattern  but  do  not  alter  the  over¬ 
all  characteristics.  If  the  ring's  size  is  doubled  and 
placed  closer  to  the  MWD  (row  3) ,  the  maximum  deflection  is 
delayed  and  somewhat  reduced,  but  not  appreciably.  There¬ 
fore,  it  is  indicated  that  the  cylinder  wall  is  still  the 
primary  resistance  medium,  and  that  the  rings  used  are  in¬ 
effective  for  controlling  initial  deformations.  A  further 
parameter  study  varying  the  ring's  cross  sectional  area, 
moment  of  inertia,  and  placement  could  better  define  the 
rings  effectiveness. 


*y  .8"  Rings 
0 . 7xl06psi 


Stringer  Stiffened  Shells  -  With  and  Without  Cutouts 

When  internal  stringers  are  added  to  the  cylinder,  the 
initial  deformation  pattern  changes.  Figure  25  depicts  a 
cylinder  with  internal  stringers  positioned  3"  on  center, 
studied  using  smeared  theory.  The  cylinder/stiffness  com¬ 
bination  is  now  able  to  resist  the  end  moments  caused  by  the 
bending  rotation  of  the  end  plane.  In  effect,  the  smeared 
stringers  have  made  the  shell  wall  thicker,  more  able  to 
act  like  a  column.  However,  the  shell  is  still  too  short 
to  allow  for  complete  flexure  analysis.  There  is  not  enough 
distance  between  the  end  constraints.  On  the  compression 
side,  the  stiffened  cylinders  follow  the  compression  zone 
radial  deflection  as  described  by  (21)  .  The  tension  side 
deforms  with  the  same  shape,  not  the  mirror  reverse  as  it 
would  be  with  total  compression  loading.  The  figure  again 
has  a  crossover  point  which  should  become  negligible  with 
an  increase  in  the  L/R  ratio.  A  study  similar  to  Figs  22 
and  23  was  not  completed,  but  it  can  be  surmised  that  the 
same  effects  would  develop;  probably  faster  since  the  shell 
is  stiffer.  Model  numbers  1,  3,  5  in  Appendix  C  were  used 
in  a  L/R  study  from  .698  to  1.361  to  contrast  results  against 
stiffened  cylinders  with  cutouts. 

Figures  26  and  27  depict  cylinders  with  1  asymmetrical 


and  2  symmetrical  cutouts  respectively.  As  can  be  seen,  the 
cutouts  affect  the  prebuckling  deformation  patterns.  In 
both  types,  the  cutouts  cause  a  commensurate  loss  in  strength 


which  lets  the  deformations  above  the  cutout  progress  as  if 
there  was  little  restraint.  The  initial  deformations  close 


to  the  top  row  follow  that  of  the  stiffened  shell  without 
cutouts,  although  the  MWD  value  is  65%  less.  However,  there 
is  no  crossover  point  anymore.  The  cutout (s)  at  the  center 
of  the  cylinder  lets  the  stresses  redistribute  to  areas  at 
the  sides  of  the  cutout.  The  bending  restraints,  which 
forces  the  crossover  to  satisfy  compatability  relations,  has 
been  altered  due  to  the  missing  section.  Hence,  the  dotted 
depiction  of  the  deflection  over  the  cutouts  assumes  almost 
a  linear  relationship.  Model  numbers  16,  18,  and  20  in  Appen¬ 
dix  C,  are  used  for  axial  compression  analysis.  The  circum¬ 
ferential  extent  is  90  degrees  but  symmetry  allows  the  full 
180°  to  be  modeled.  However,  this  lets  us  only  examine  the 
case  of  two  cutouts  under  axial  compression.  For  bending, 
model  numbers  17,  19  and  21  A  4  B,  are  used.  Minor  input 
card  changes  in  the  subroutine  wall  allow  us  to  consider 
either  one  (B)  or  two  cutouts  (A)  respectively. 

The  stiffened  cylinders  were  loaded  with  the  same 
bending  rotation  angle  as  was  used  for  the  unstiffened 
analyses.  The  stress  intensity  at  the  90°  column  (maximum 
compression  column)  of  row  1  at  buckling  was  contrasted 
against  the  same  value  under  axial  compression.  Figure  28 
depicts  the  results.  The  general  movement  for  the  entire 
graph  is  progressively  higher  with  increasing  L/R.  This  im¬ 
plies  that  the  buckling  loads  increase  with  length.  How¬ 
ever,  my  analysis  runs  were  for  limited  changes  of  length 
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of  short  cylinders.  Nelson  (21)  plotted  the  same  pro¬ 
gressive  movement  for  SS4  boundary  conditions  but  went  on 
to  show  that  for  a  greater  L/R  ratio  equal  to  1.78  the 
graph  drops  down  again  to  below  the  values  for  L/R  =  .698. 
The  pertinent  mechanism  at  work  is  the  interaction  of  the 
SS4  boundary  condition  with  the  short  shell.  Nelson  shows 
that  for  other  simple  support  (SSI,  SS3)  boundary  re¬ 
straints,  the  plot  remains  close  to  horizontal.  As  was 
brought  out  earlier  in  the  shell  model  section,  we  are  re¬ 
quired  to  use  the  SS4  boundary  support  for  analysis,  due 
to  STAGSC  handling  of  displacement  loaded  buckling  analy¬ 
ses.  Pure  bending  rotation  is  a  form  of  displacement 
loading.  Therefore,  this  boundary  condition  will  affect 
the  comparison  between  axial  and  rotation  loading  analyses. 

For  stiffened  cylinders  without  cutouts.  Fig  28  shows 
that  the  bending  critical  values  are  above  the  axial  com¬ 
pression  loaded  values  and  tend  to  diverge  slightly  as  the 
L/R  ratio  increases.  An  explanation  for  this  behavior, 
involves  crossover  point  analysis  under  prebuckling  dis¬ 
placement  rotations.  As  was  pointed  out  earlier,  an  inter¬ 
nally  stringer  stiffened  cylinder  deforms  in  the  shape  of 
Fig  25.  When  axial  compression  is  compared  with  the  shape 
for  the  same  L/R,  the  axial  displacement  pattern  at  90° 
follows  the  same  concave/convex  arrangement.  However,  the 
crossover  point  is  closer  to  the  top  edge  for  the  axially 
compressed  cylinder.  The  variation  is  probably  due  to  the 
effect  of  stringer  eccentricity  aiding  the  initial  defor- 
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formation  under  bending  and  resisting  the  return  mechan¬ 
ism.  The  difference  in  crossover  points  causes  a  commen¬ 
surate  change  in  the  effective  length  (distance  between 
the  cylinder  top  half  and  bottom  half  crossover  points) 
between  the  axial  and  bending  loaded  models.  Since  the 
bending  effective  lengths  are  smaller  in  each  case,  their 
buckling  loads  are  slightly  higher.  This  effective  length 
difference  increases  slightly  with  an  L/R  increase  which 
causes  the  divergence.  These  effects  are  associated  with 
short  shells  only;  as  soon  as  "beam  type"  flexure  predom¬ 
inates,  the  stress  analysis  models  change.  However,  for 
these  cases  the  maximum  difference  was  less  than  13%. 
Therefore  bending  buckling  stresses  can  be  modeled  by 
axial  compression  analysis  to  within  that  percentage  tol¬ 
erance  for  the  given  L/R  ratios. 

The  picture  changes  when  cutouts  are  introduced.  Now 
there  is  less  material  under  both  loading  schemes  and  the 
buckling  intensities  are  lower  for  all  L/R  ratios  compared 
with  no  cutout  cylinders.  Further,  the  maximum  bending  loads 
have  less  material  to  average  over  and  consequently  the  cyl¬ 
inder  buckles  earlier  when  compared  against  axial  compression 
analysis.  Although  the  SS4  boundary  condition  is  still 
present,  the  cutout  condition  is  now  dominant.  But,  as  the 
shell  length  increases,  the  hole  (size  held  constant)  effect 
diminishes  in  importance  and  the  critical  axial  vs  bending 
intensity  values  converge  between  themselves  and  the  no  cut¬ 
out  values. 
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Two  other  points  can  be  made  from  Fig  28.  First, 
there  is  negligible  difference  between  the  one  cutout  and 
two  cutout  models  under  bending  loading.  The  critical 
values  are  within  .1%  of  each  and  therefore  plotted  to¬ 
gether.  This  effect  is  also  clearly  demonstrated  on  Fig 
32  where  the  one  cutout  and  two  cutout  stress  intensity 
flow  lines  lie  right  beside  each  other.  The  reason  for 
the  unimportance  of  the  other  cutout  is  that  it  lies  in 
the  tension  zone,  which  does  not  contribute  to  the  maximum 
compression  fiber  analysis.  Second,  a  "finer"mesh  (15  x 
121)  is  plotted  on  the  figure  for  axial  compression  analy¬ 
sis  of  a  stiffened  cylinder  without  cutouts.  For  unstiff¬ 
ened  cylinders,  increasing  the  mesh  from  15  x  61  to  15  by 
121  resulted  in  a  buckling  load  increase  of  19%.  For 
stiffened  cylinders,  this  same  mesh  refinement  only  brings 
out  an  increase  of  5.8%.  This  indicates  that  for  stiffened 
analysis  a  moderate  spacing  models  the  cylinder  very  well. 
This  was  the  justification  for  using  the  15  x  61  mesh 
arrangement  as  the  primary  strategy. 

Figure  29  depicts  the  critical  bending  moment  values 
for  the  different  cutout  and  no  cutout  cylinder  models  as 
the  length  is  changed.  The  algorithm  used  to  calculate 
the  moment  is  given  in  Appendix  B.  Since  the  cylinder  is 
loaded  by  specifying  a  rotation  angle,  the  applied  moment 
is  found  as  a  by  product  from  the  resulting  Rvalues.  The 
SS4  boundary  condition  again  causes  the  rise  in  the  criti¬ 
cal  moment  value  as  the  L/R  ratio  increases  for  these 
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short  cylinders.  What  is  interesting  is  the  closeness  of 
the  different  critical  moment  values  for  the  shell  cylin¬ 
der  models  in  question  for  each  L/R  ratio.  This  indicates 
that  we  are  able  to  predict  critical  moment  cutout  values 
using  no  cutout  model  analyses.  The  maximum  difference 
lies  within  9%  for  the  L/R  range  studied.  There  is  some 
difference  between  the  one  cutout  and  two  cutout  values 
which  arises  from  using  the  top  row  values  for  the  total 
model  to  calculate  the  total  moment.  However,  this  differ¬ 
ence  is  always  less  than  9 X. 

Plots  of  the  eigenvectors  for  the  no  cutout  and  sym¬ 
metrical  cutout  stiffened  cylinder  models  are  given  in 
figures  30  and  31  respectively.  These  figures  are  included 
to  show  that  the  circumferential  buckling  mode  shape  under 
end  rotation  loading  doesn't  change  with  the  addition  of  a 
cutout.  The  axial  compression  loading  mode  shapes  change 
but  only  in  the  relative  magnitude  of  displacements.  For 
all  cases,  the  wave  lengths  are  the  same,  giving  approxi¬ 
mately  12  half  sine  waves  per  180°  compression.  The  axial 
wave  form  for  both  the  cutout  and  non  cutout  stiffened 
models  is  one  %  sine  wave,  with  the  maximum  normalized  de¬ 
flection  occuring  at  the  center  of  the  shell.  This  is 
contrasted  with  unstiffened  cylinders  which  have  many  half 
sine  waves  in  the  axial  direction. 

An  attempt  was  made  to  correlate  the  cutout  model 
loading  under  axial  compression  to  that  under  pure  end 
bending  by  means  of  a  ratio  depicted  in  Fig  33.  The 


approach  was  motivated  by  the  convergence  of  the  bending 
vs  axial  values  in  Fig  28  and  the  equal  number  of  sine 
waves  over  the  cutout  for  the  different  loading  schemes. 
Figure  31.  It  was  felt  that  the  total  load  above  the  cut¬ 
out  could  be  used  as  a  predictive  technique.  The  total 
load  above  the  cutout  for  both  loading  schemes  was  arrived 
at  by:  1.  Multiplying  each  top  row  prebuckling  Nx  value 
above  the  cutout  by  the  distance  it  works  over.  2.  Multi¬ 
plying  each  of  the  above  values  by  the  eigenvalue  at  buck¬ 
ling.  3.  Summing  the  values.  These  total  loads  were  then 
used  in  the  ratio  in  Fig  33  and  graphed.  The  result  was  a 
plot  that  decreased  as  the  length  increased  which  was  ex¬ 
pected  due  to  the  convergence  of  Fig  28  points.  Expected 
values  for  bending  buckling  loads  above  a  cutout  can  be  ob¬ 
tained  using  these  ratio  values,  especially  for  L/R  >  1. 


W  Deflection  V  Deflection  W  Deflection 

Figure  25.  Prebuckling  Deflections  for  Stiffened 

Cylinder 


W  Deflection  V  Deflection  W  Deflection 

Figure  26.  Prebuckling  Deflections  for  Stringer 
Stiffened  Cylinder  With  1  Cutout 
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Figure  29.  Buckling  Moment  vs  Length  For 
Stiffened  Cylinders  With  an3 
Without  Cutouts  Under  Pure  End 
Rotation 


50|tr*  >  01  ft  M 


Compression  Half 


31.  Eigenvectors  For  Axial  Compression  and 
Pure  Bending  Loading  of  a  Stiffened 
Cylinder  With  2  Symmetrical  Cutouts 


—  Cutout  - #( 


Figure  32.  Stress  Intensity  Flow  For  Axial  Compression 
and  Pure  Bending  Loading  of  a  Stiffened 
Cylinder  at  Buckling 


Where  the  Critical  Load  Cutout  Factor  (CLCF)  is 


Axial  Load  _  Bending  Load 
Above  Cutout  Above  Cutout 


Axial  Lead  Above  Cutout 


With  the  Loads  Calculated  Using  Nx  Values  Along 
Row  One  Above  the  Cutout  at  Buckling. 


Figure  33.  Calculated  Loading  Ratio  Above 
Cutout  vs  Change  in  Length  For 
A  Stiffened  Cylinder 
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V.  Conclusions 


From  the  analysis  undertaken  in  this  study,  the  follow¬ 
ing  conclusions  can  be  made. 

(1)  The  Skogh-Brogan  routine  for  applying  pure  bending 
rotations  to  a  cylinder  using  end  boundary  kinematic 
constraints  does  keep  the  edge  circular  and  planer 
under  loading. 

(2)  When  the  pure  bending  loads  are  applied  to  short  un¬ 
stiffened  cylinders  (L/R  <  10)  the  structure  assumes 
deformation  patterns  different  from  long  cylinder 
bending  deflections.  The  material  on  the  tension 
zone  (-90°)  deforms  inwards;  whereas,  the  material 
on  the  compression  zone  (+90°)  bulges  outwards. 

These  radial  displacements  are  of  equal  magnitude 
but  opposite  in  sign.  The  effect  diminishes  cir¬ 
cumferentially  to  zero  at  the  0  and  180°  points. 

The  longitudinal  extent  of  this  deformation  pattern 
changes  both  by  degree  and  direction  with  cylinder 
length. 

(3)  The  inplane  V  deflections  at  the  middle  columns  (0° 
and  180°)  follow  a  "beam  type"  deflection  path 
reaching  a  maximum  at  the  center  mid  plane  of  che 
cylinder . 

(4)  As  the  cylinder  length  increases  the  effects  dimin¬ 
ish  until  an  L/R  =  20  where  there  is  no  discernible 
deformation  for  the  mesh  spacing  used.  The  de- 
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flections  at  this  length  follow  "beam  type" 
flexure. 

(5)  Discrete  ring  stiffening  at  the  ends  of  the  cylin¬ 
der  did  not  appreciably  alter  the  overall  charac¬ 
teristics  of  the  deformations;  at  least  for  the 
chosen  stringer  dimensions  and  positioning  from  the 
loaded  edge. 

(6)  When  stringers  were  added,  the  deformation  pattern 

was  altered.  The  cylinder/stringer  combination  is 
stiff  enough  at  the  boundary  to  follow  beam  flexure 
behavior.  However,  with  shorter  lengths  L/R  <  10 

there  isn't  enough  material  between  the  edges  to 
allow  for  this  flexure.  Therefore,  the  cylinder's 
radial  displacements  follow  a  reversal  of  direction 
near  the  center  mid  section.  However,  the  inplane 
V  deflections  at  the  middle  column  continues  to 
follow  the  typical  flexure  pattern. 

(7)  When  cutouts  are  added  to  the  stiffened  cylinder  un¬ 
der  bending  loading  at  the  maximum  compression  and 
tension  zones,  the  deformation  pattern  initially 
follows  that  for  no  cutout  stiffened  cylinders.  As 
the  cutout  is  approached,  the  cylinder  loses  its 
rigidity  and  the  deformations  become  almost  linear. 

(8)  For  buckling  analysis,  the  tension  zone  has  very 
little  effect  on  the  critical  stresses  of  the  maximum 
compression  zone  fibers.  Adding  a  cutout  in  the 


tension  zone  does  not  alter  the  failure  values. 


(9)  The  buckling  stress  intensity  values  for  stiff¬ 
ened  cylinders  without  cutouts  under  axial  vs 
bending  loading  lie  within  13%  of  each  other  for 
the  L/R  range  of  .698  to  1.361.  Axial  analysis 
can  be  used  therefore  to  predict  moment  load  fail¬ 
ures  to  at  least  this  percentage  accuracy. 

(10)  Bending  vs  axial  analysis  for  stiffened  cylinders 
with  cutouts  is  less  precise.  The  trend  shows  a 
convergence  of  critical  values  to  within  9%  as  the 
length  increases  to  L/R  =  1.361.  For  L/R  values  >  1, 
a  load  ratio  graph  along  with  axial  data  can  be 
used  to  model  bending  loads  above  the  cutout. 
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Appendix  A 


Force  and  Moment  Intensities .  Constitutive 
Equations ,  and  Kinematic  Relations 

The  internal  forces  and  moments  for  a  shell  are 
usually  expressed  in  terms  of  forces  and  moments  per  unit 
distance  along  the  edge  of  the  shell  element.  These  force 
and  moment  intensities,  pictured  on  Fig  1,  are  related  to 
the  stresses  as  shown  on  the  next  page.  The  symbols  a*  and 
Tjj0  denote  stress  components  at  any  point  through  the  shell 
wall  thickness. 

The  following  page  contains  the  kinematic  middle- 
surface  relations  postulated  by  Donnell  and  the  constitutive 
equations  for  thin-walled  isotropic  elastic  cylinders.  The 
form  of  these  equations  is  taken  from  (14) . 


Appendix  B 


STAG SC  PROGRAM  LISTINGS 

The  following  listings  show  the  different  aspects 
of  the  STAGSC  computer  code  I  used.  They  do  not  come  from 
the  same  program. 

B1  is  a  listing  of  the  input  data  cards  needed  to  run 
a  program.  This  run  was  for  a  cylinder  of  total  length 
40  inches  with  a  mesh  arrangement  of  16  by  71  over  the 
quarter  shell  segment  model.  It  has  options  included 
for  bifurcation  analysis  of  a  dummy  branch  moment  loaded 
cylinder  with  a  wall  subroutine  defined  surface  with  12 
inch  cutouts. 

B2  details  the  user  written  subroutine  UCONST  for  a 
cylinder  segment  having  61  columns. 

B3  shows  the  user  written  subroutine  WALL  with  no 
cutouts . 

B4  shows  the  user  written  subroutine  WALL  with  two 
cutouts . 

B5  details  the  subprogram  MOMENT  which  calculates 
the  total  applied  moment  from  the  Nx  values  along  the 
top  row  of  the  cylinder. 


INPUT  DATA  CARDS 
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SUBROUTINE  UCONST 

PURE  MOMENT  -  EDGE  STAYS  CIRCULAR  AND  IN  CNE  PLANE 
INPUT  IS  U  DISPLACEMENT  AT  ROW1/CGL1  ON  DUMMY  BRANCH  2 

OX  HENS  I  ON  IBRNCH(S),  IX  (5),  IY(S),  ID<5>,  CC<5>,  HA(Si) 


OATA  T-U/-90., 

*8  7  •  y  * 8  "4 

. y  “81. y  “ 76  *  y 

-71 .  , 

“72., 

*69. 

y  “66. y 

“63.  y  -60#  y 

“E 7.  y 

• y  “51. y  “ns . y 

“4£  .  . 

“H2.y 

-39. 

y  “35. y 

“33.  y  "3  O’,  y 

“27 •  y  -24 

• y  “21. y  “18 • y 

-If 

-12.  y 

-9., 

“6  .  y  “3  .  y 

O.Uy  3. y  & • y 

* 

m 

CM 

-H 

* 

• 

0> 

15  •  y  18  ■  t  2 1.  y 

2A.y 

27., 

3)  .  y 

33., 

C  36ii  3S«t  42.,  45.  y  s  fiti  Sid  64.  y  57»>  bO.y  63«,  bS.y  b9»  , 

C  72»,  75t  i  7 8*  •  81.  y  64.,  87.,  90*/ 

E  *  IO.4E+6 
IX(i)  «  1 
IX  (3)  s  1 

oo  io  r=i,6i 

TH  =  THA(I)*3«lt  159  256/1 60* 

IX (2)  *  1 
IBRNCH(l)  =  1 
IBRNCHI 2)  *  1 
1BRNCH ( 3 )  *  1 

THIS  SECTION  PREVENTS  SINGULARITIES 
IF (I«E2« 31)  GO  TO  5 
XF(I.EO.l)  GO  TO  2 
ZFiI*E2.Sl>  GO  TO  2 

EQUATION  1 
IY41)  *  I 
10(1)  *  2 

CC(1)  «  -SIN(TH)  *  E 
IY (2)  *  I 
10(2)  *  3 

CC (2)  *  COS (TH)  *  E 

CALL  SONSTR  (2,IBRNCH,IX ,IY,IO,CC) 

2  CO NT  IN JE 


B2 
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UCONST  cont. 


EQUATION  2 
IY(1)  *  31 
10(1)  *  2 
CCC1)  *  ♦£ 

IYC2)  *  I 
10(2)  =  2 

CCC2)  *  -C3S(TH)*E 
IY (3)  *  I 
10(3)  «  3 

CC(3)  *  -SIN(TH)*E 

CALL  20NSTR  (3.IBRNCH ,IX » IY. I0,C5) 

EQUATION  3 
IY (1)  *  31  • 

10(1)  =  1 
CC(i)  =  E 
IBRNCHC 2)  »  2 
IX  (2)  *  1 
IY  (2)  *  1 
10(2)  *  1 
CC(2)  *  SIN (TH) *E 
IY (3)  *  I 
10(3)  *  1 
CC(3)  *  -E 

CALL  3DSSTR  (3tIBRNCMfIX »1YiID»C3> 
5  CONTINJE 
10  CONTINUE 
RETURN 
ENO 
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WALL  -  NO  CUTOUTS 


SUBROUTINE  WALL(IBRNCH,X,Y,Z,ICF9,ISTFP,IPRW,RH0A,C:S) 
COMMON/lAYDl/  TL  (20>,tX3  (20)  .EY3C20)  ,U2i(20>  ,G3  ( 20  »  ZET  3  <20 ) . 
1RH03(23) » LAYS. LSTRS 

COMMON/SNEAR/NST I» ITYP<5 ) »IDIR( 5) » SPACE  (5) ,ANGE(5) 

C 

IF  (IBRNCH.NE.l)  RETURN 
LAYS  *  1 
1CFB  *  3 
IPRH  =  0 
C  PANEL  DATA 

TL (1)  =  0.1 
EX 3(1)  *  lOAOCOOC. 

EY 3 ( 1 )  *  1U40000 (i. 

U21(  1)  =  .33 
03(1)  s  3750000. 

ZET3(1 )  *  0.0 
RH03(1)  *  0.1 
C  SMEARtO  DATA 
ISTFF  *  1 
NSTI  *  L 
ITYP(l)  «  1 
IOIR(l)  *  1 
SPACE(l)  =  3.0 
RETURN 
ENO 
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WALL  -  TWO  CUTOUTS 


SUBROUTINE  WALL( IBRNCHyX  »Y»Z»ICF3* ISTFF  ,IPRH,RHOA,C: 2) 
COMMON/. A YOl/  TL (20) » EX3 (20) »EY3<20)  ,U21(20)  ,G3 C 20) » ZET 3 <2C ) , 
1RH03(20) ,LAYS,LSTRS 

COMMON/ 5 M EAR/NS T I »ITYP(fe  }  » 1GIR(5)  •  SPACE  IB)  ,ANGE(5) 

IF  <IBRNCH<NE< 1)  RETURN 
LAVS  a  1 
ICFB  *  3 
IPRH  *  9 

CUTOUT  LOCATION  TESTS  FOR  12BY12  INCH  CJT 
BOTTOM  LE?T 

IF  (X  ,3T.  14.03  < AND <  Y  <LT<  -84.3->  GO  TO  13 
BOTTOM  RIGHT 

IF  (X  .ST#  14.  Qo  .AND.  Y  <GT  «  b+.Jl)  GO  TO  10 
PANEL  DATA 
TL  (1)  =  0.1 
EX3( 1)  s  10400000. 

£Y3(i)  3  104000CC. 

U2 1 ( 1 )  s  .33 
G3 (1)  3  375C000. 

ZET3I1)  3  0.0 
RH03(1)  3  o,i 
SMEAREO  DATA 
ISTFF  s  1 
NSTI  s  1 
ITYP<11  3  1 
IOIR(l)  3  1 
SPACECl)  3  3.0 
GO  TO  30 
10  CONTINJE 
CUTOUT  DATA 
ISTFF  3  0 
TL  (1)  3  0,1 
£X3(1)  3  o,Q 
EY3C1)  =  0.0 
U21(l)  3  G •  0 
G3 (1)  3  0.0 
ZET3<1)  b  q,o 
RH03(1)  3  0,0 
30  CONTINUE 
RETURN 
END 
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PR0GR4  H  MOMENT  « NPUT,  OUT  PUT « TAP£3=IMPUT  ,TA PE4«0UTPUr ) 

2  REAO  (3,1)  A 

1  FORMAT  ( 25X  * A2) 

IF  (A.NE.2HNX)  GO  TO  2 
SMOMENT  =  0,0 
00  10  I  =  1,91 
REAO  (3,3)  A,X 
WRITE  (4,3)  A ,X 

3  FORMAT  (8X,F8.4,- X,E12.6) 

AA  =  A*  3.14159256/180 •& 

IF  (  I  ,EQ.  1  )  GO  TO  21 

IF  (  I  *  GT •  1  .AND.  I  .IT,  31)  GO  TO  30 

IF  (  I  • EO.  31  ) GO  TO  41 

IF  (  I  ,EQ.  91  )  GO  TO  9 C 

S  *  X  *  1.5  *  SIN(AA) 

GO  TO  S3 

20  S  *  X  *  1.6  *  SIN(AA) 

GO  TO  3) 

30  S  a  X  *  3.0  *  SIN(AA) 

GO  TO  30 

40  S  *  X  *  2.25*  SIH(AA) 

GO  TO  S3 

50  S  *  X  *■  .75  *  SIN(AA) 

60  SMOMENT  *  SMOMENT  ♦  S 
10  CONTINJE 

XX  *  2,3  *  57.2958  *  SMOMENT 
WRITE  (4,4)  XX 

4  FORMAT  ( E12.6) 

STOP 

ENO 
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APPENDIX  C 


FINITE  DIFFERENCE  MESH  ARRANGEMENTS 

The  different  mesh  arrangements  are  by  model  number. 

Models  1  to  15  and  2A  pertain  to  unstiffened  and  stiffened  shells 
without  cutouts.  Models  16  to  21A,  B  are  used  for  the  cutout 
analyses  only. 

Model  3,  the  15  by  61  mesh  arrangement  is  the  basic  model 
around  which  the  other  non  cutout  models  are  compared.  It  gives 
a  3"  by  3"  square  grid  spacing  for  most  of  the  cylinder  segment 
with  a  1"  by  3"  ratio  for  the  top  few  rows.  For  the  -.iffened 
analyses,  it  gives  results  with  a  circumferential  buckle  wave 
pattern  of  4  nodes  per  half  sine  wave,  without  having  a  grid  line 
between  each  stiffener  (economical).  For  unstiffened  analysis  it 
gives  consistent  although  conservative  buckling  load  results, 
approximately  78%  of  the  classical  value.  Since  the  unstiffened 
analyses  were  mainly  shape  studies,  these  results  are  close 
enough.  Models  8  to  15  were  used  for  moment  end  deflection  analysis 
only  and  not  used  for  bifurcation  studies. 

Cutout  models  16,  18,  20  are  1/8  shell  Nelson  models  used  for 
stiffened  axial  compression  analysis.  Models  17,  19  and  21  are 
used  for  the  bending  studies.  The  A  and  B  suffixes  refer  to  two  cuts 
and  one  cut  respectively.  The  dashed  cutout  is  the  one  not  there 
for  one  cut  analysis. 


96 


99 


ADDITIONAL  MODELS  FOR  SHELLS  WITHOUT  CUTOUTS 


MODEL 

SHELL  WIDTH 

Y  DIRECTION 

MESH 

NODES 

2A 

90" 

30  at  3" 

465 

Note : 

x  direction  mesh  remains 

the  same  for 

shell 

length  of  60" 

MODEL 

SHELL  LENGTH 

X  DIRECTION 

MESH 

NODES 

8 

120" 

6  at  1" 

9  at  6" 

496 

9 

180" 

5  at  1" 

17  at  5" 

713 

10 

240" 

5  at  1" 

23  at  5" 

899 

11 

300" 

6  at  1" 

24  at  6" 

961 

12 

360" 

5  at  1" 

25  at  7" 

961 

13 

480" 

6  at  1" 

26  at  9" 

1023 

14 

600" 

3  at  1" 

27  at  11" 

961 

15 

1143" 

3  at  1" 

15  at  6" 
48  at  12" 

1829 

Note:  Y  direction  mesh  remains  the  same  30  at  6" 


6" 

6  at  1" 
l 


*-15"  - 

LOatl.5 


r 


_  75"  _ 

25  at  3" 


M 


MODEL  20  23  by  36  (1  Cut) 


6"  ►”15"  — * 

6  a|t  1M  10atl.5 


150" 


50  at  3" 


h-  15"-* 
10atl.5 
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